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Improvements in testing and modeling of nonlinear and unsteady aerodynamic effects for 
flight dynamics predictions of vehicle performance is critical to enable the design and 
implementation of new, innovative vehicle concepts. Any configuration which exhibits 
significant flow separation, nonlinear aerodynamics, control interactions or attempts 
maneuvering through one or more conditions such as these is, at present, a challenge to test, 
model or predict flight dynamic responses prior to flight. Even in flight test experiments, 
adequate models are not available to study and characterize the complex nonlinear and 
time -dependent flow effects occurring during portions of the maneuvering envelope. 
Traditionally, airplane designs have been conducted to avoid these areas of the flight 
envelope. Better understanding and characterization of these flight regimes may not only 
reduce risk and cost of flight test development programs, but also may pave the way for 
exploitation of those characteristics that increase airplane capabilities. One of the hurdles is 
that the nonlinear/unsteady effects appear to be configuration dependent. This paper 
compares some of the dynamic aerodynamic stability characteristics of two very different 
configurations - representative of a fighter and a transport airplane - during dynamic body- 
axis roll wind tunnel tests. The fighter model shows significant effects of oscillation 
frequency which are not as apparent for the transport configuration. 
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I. Introduction 

Improved prediction of airplane aerodynamic characteristics is required to improve design performance and 
safety, while at the same time reducing development and certification costs. Missed predictions of aerodynamic 
effects on flight dynamics and handling qualities of airplanes continue to give rise to large schedule impacts and 
costs when discovered in flight test. Aviation history is replete with such examples during the development of 
nearly every tactical fighter airplane 1 . Usually the large discrepancies between predicted and actual characteristics 
discovered in flight occur in non-linear or maneuvering flight regimes, such as high angles of attack, transonic 
maneuvering, or high angular rate maneuvering . For commercial airplanes, it is becoming more important to predict 
flying qualities in off-nominal flight conditions as well. Improved predictions may enable better departure resistance 
and recovery strategies, and may provide for higher fidelity pilot training in unusual flight conditions. Enhanced 
aerodynamic models for extreme envelope conditions will also improve analysis of airplane accidents and incidents, 
thereby providing a basis for improving safety of the airplane fleet. 

One of the primary areas in which aerodynamic models are known to be deficient is in their response to high 
dynamic rates. Traditionally, dynamic rate effect models have consisted of coefficients that are linear functions of 
rate as obtained, classically in derivative form, from forced oscillation testing, and coefficients that are nonlinear 
functions of wind-axes axes roll rate as obtained from rotary balance testing. The resulting aerodynamic model 
structure computes the forces and moments based on instantaneous values of angular rates and the flow incidence 
angles, a and p. Forced oscillation results for many configurations often show large variation in the linear damping 
derivatives with oscillation amplitude and frequency, especially near stall angles of attack. These variations are 
indicators of either nonlinear rate effects or unsteady aerodynamic effects. Over the past few years as data handling 
and storage capabilities have improved, the time history data of complete forced oscillation cycles have been used 
for analysis, providing far more information than just the traditional linear derivatives. This research has led to 
improved understanding of the nonlinear and time-dependent effects present during unsteady airplane motion. In 
turn this has led to the development of higher fidelity modeling capabilities of airplane flight dynamics 2 3 " 6 . 

Unfortunately, there appears to be a large variation in aerodynamic response depending on the airplane 
configuration tested, which adds to the difficulty of developing a proper physics-based model to describe airplane 
dynamic characteristics. This paper will compare the aerodynamic response to roll oscillations from low-speed 
forced oscillation wind tunnel tests of two very different configurations: a high-aspect ratio configuration 
representing a modern transport airplane, and a low-aspect-ratio configuration representing a modern fighter 
airplane prototype. Both sets of data were measured on the same testing apparatus and in the same wind tunnel using 
the same data reduction algorithms. Although data were obtained during oscillations in all three axes (pitch, roll, 
yaw), only roll axis data will be presented here. Additionally, many configuration variations were tested including 
control surface deflections, component build-ups, etc.; however, only one configuration for each airplane type will 
be presented. 


II. Wind Tunnel Tests 

Sketches of the two configurations tested are shown in Figs. 1 and 2. The fighter model (Fig. 1) features a 
blended wing/body with an aspect ratio of 2.1, twin vertical tails and no horizontal tails. The transport model (Fig. 
2) is a twin-engine configuration and has a slender fuselage, slender wings with an aspect ratio of 7.9 and a 
conventional horizontal and vertical tail arrangement. Key configuration details are listed in Table 1. 
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Table 1. - Geometry Comparison of Models 


Parameter 

Transport Model 

Fighter Model 

Wing Span, ft 

6.85 

5.40 

Reference Wing Area, ft" 

5.90 

13.95 

Aspect Ratio 

7.9 

2.1 

Mean Aerodynamic Chord, ft 

0.92 

3.3 


The fighter model is shown mounted in the NASA Langley Research Center’s 14x22-Foot Subsonic Wind Tunnel in 
Fig. 3. The transport model is shown mounted in the same wind tunnel in Fig. 4. Data were taken over a wide range 
of conditions, including angles of attack from ±90°, oscillation amplitudes from ±5° to ±30°, and oscillation 
frequencies from 0.23 to 1.10 Hz, resulting in reduced frequencies ( k = a>b/2V) of k = 0.054 to k = 0.455. Table 2 
lists the range of data obtained that will be presented in this paper. 


Table 2. - Data Envelopes for Comparisons 


Parameter 

Transport Model 

Fighter Model 

Amplitudes of oscillation 

0, ±5°, ±10°, ±20°, ±30° 

0, ±5°, ±15°, ±30° 

Oscillation Frequency, Hz 

0.23-0.92 

0.3-1.10 

Reduced Frequency, k 

0.054-0.216 

0.055 - 0.455 

Dynamic Pressure, psf 

10 

2, 5, 10 

Angle of attack 

0° - 90° 

-5° - 75° 

1 p 1 maximums during oscillations 

0.018-0.113 

0.005 - 0.238 


The aerodynamic components of the dynamic data were obtained by first oscillating the model at set 
frequencies and amplitudes with the wind off to obtain gravity and inertial tare data and then by repeating the 
process with the wind on and subtracting the tare data. During both the wind-on and wind-off runs, data were 
obtained over 40 oscillation cycles, sampled at 200 Hz with 100 Hz anti-aliasing filters. The resulting data were then 
further filtered with a 4-Hz lowpass filter to remove unwanted frequency content. The dynamic data could then be 
ensemble averaged to represent the variation in aerodynamic forces and moments over a single oscillation cycle. 
Concurrently, the traditional calculations of damping derivatives using the in-phase and out-of-phase components of 
the balance outputs were calculated at the conclusion of each set of oscillations cycles. The computed damping 
derivatives actually combine roll-rate-damping and sideslip-rate effects in one single parameter, 
Cj p = C tpi + C { sin(a) . In addition to the dynamic data, complementary static wind tunnel data were obtained on 

the same test rig for comparison. 

An indication of the “linear” characterization of forced oscillation data can be gained by inspection of the 
moment coefficient vs. oscillation angle. A conventional modeling approach C t = C) (5 + Q p yields an inclined 

oval as depicted in Fig. 5. The inclination of the oval represents the static stability, or restoring moment due to roll 
angle (or sideslip) and the vertical diameter of the oval represents the damping effect, or restoring moment due to 
angular rate. Deviations from this oval shape indicate non-linear effects that cannot be modeled using the summation 
of linear, first-order derivatives. In order to capture these non-linear effects, higher order functions or table lookups 
are often required. 


III. Fighter Model Data Results 

Static data (Fig. 6) show maximum lift for the fighter model occurs at approximately a = 34°. with an earlier 
decrease in the lift curve slope noticed at a = 24°. Static rolling moment characteristics (Fig. 7) show a large 
asymmetry which also occurs near a = 24°. This combination of characteristics would certainly warrant further 
investigation as to the mechanisms and dynamic characteristics if flight near this angle of attack was anticipated. 
Figure 8 shows the classical results from forced oscillation testing. For this example, forced oscillation tests with an 
oscillation amplitude of ±30° were conducted at a number of reduced frequencies, k. The test data show that at low 
a’s, roll damping (Cjp) is relatively invariant with oscillation frequency suggesting that the roll damping parameter 
is also apparently invariant across the range of roll rates tested, since at a given amplitude, varying frequency also 
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varies the maximum roll rates generated. At higher angles of attack, the data is not as well behaved as evidenced by 
large variations of the roll damping parameter with frequency from approximately a = 24° to a = 40°. This 
characteristic behavior has been observed in tests for decades, and until recently, has only been addressed by 
selecting model test frequencies near full-scale airplane response frequencies of interest (i.e. dutch-roll frequency 
for example), or by parametrically varying roll damping across a wide range of potential values to ensure sufficient 
robustness in the flight control system design. 

More recently, testing and analysis has been conducted in an attempt to reduce the uncertainty in the dynamic 
characteristics measured in forced oscillation tests. An alternative data reduction approach, called the single-point 
method, uses data obtained from time history recordings of the forced oscillation cycles. This method uses 
instantaneous values of the aerodynamic forces and moments measured as the model moves through the midpoint of 
the sinusoidal rolling motions where the oscillation amplitude is zero, p = 0°, rotation rate is a maximum, and there 
is no rotational acceleration. The value of the aerodynamic forces and moments at these points (one at a positive 
rate, and one at a negative rate) can then either be divided by the non-dimensional rate at each of the points to obtain 
a derivative, or can be plotted versus rate to develop a nonlinear dynamic model that is a function of rate The 
variation of nondimensional roll rate can be accomplished by changing either the frequency or the amplitude of 
oscillation. To maintain constant reduced frequencies for rate variations, both oscillation frequencies and amplitudes 
must be varied. The test technique is limited to variations of two of the critical parameters ( p ,k, amplitude) in order 
to keep the other constant. Therefore, amplitude effects - due to transitions from one flow-field state to another 
during the oscillation - and frequency effects, which are an indication of unsteady aerodynamic phenomena, may be 
difficult to isolate. 

Figures 9-11 show the rolling moment characteristics as the model is oscillated. The cases were selected such 
that each data set was measured at the same nondimensional roll rate ( p = 0.04) as the model rolled through the 
midpoint of the oscillation, at <f) = 0°. For a linear, symmetric model of the roll damping to be valid, data from each 
amplitude of oscillation should coincide at the <f> = 0° point, since the states (a, p, p ) are identical for each data set. 
At a = 10° (Fig. 9), well before stall where the aerodynamic responses are primarily linear, the data all coalesce at 
(J) = 0° regardless of amplitude of oscillation. As a is increased to a = 30° (Fig. 10), near the stall angle of attack, 
the data loops do not coalesce at <f) = 0°, and show that as the oscillation amplitude increases, there is an apparent 
increase in roll damping. At post-stall a’s (Fig. 11), again the data do not coalesce at (f) = 0°, however the smaller 
amplitude data show higher apparent damping for this condition. The observed nonlinearity of the results at the stall 
and post-stall a conditions indicate a higher order model is needed to describe the data. 

Figures 12-14 show single-point data versus nondimensional roll rate. The plots include data obtained from a 
range of oscillation amplitudes and frequencies, at several wind tunnel velocities. Typically, three points along the 
pb/2V- axis are plotted for each reduced frequency, each representing different oscillation amplitudes, to show the 
variation with rate at a given reduced frequency. Previous comparisons of the dynamic roll data have shown tunnel 
velocity effects to be insignificant, with reduced frequency and rate being the predominant effects. The data in Fig. 
12 show that at a = 10°, the rolling moment variation with roll rate is slightly nonlinear, but apparently no 
significant reduced frequency effects are evident. At a = 30° (Fig. 13), there is a clear indication of the variation in 
the dynamic stability characteristics with reduced frequency. These data correlate well to the data in Fig. 8, which 
also show variations in roll damping with reduced frequency. At very high a’s, the data show pronounced 
nonlinearities in both reduced-frequency and roll-rate effects as seen in Fig. 14 at a = 50° . The data in Fig. 14 
show roll damping is unstable at low rates, and stable at high rates. This roll damping characteristic can give rise to 
wing rock type motions if the static stability ( C/ ) provides the required “spring action” to the dynamic system. 

IV. Transport Model Data Results 

The transport model was tested in a flaps-up cruise configuration. Static data (Fig. 15) show initial flow- 
separation wing stall occurring around a = 10°, followed by a second post-stall increase in planform lift up to a 
maximum at approximately a = 40°. The flow separation-induced stall manifests itself as a discontinuity in the lift 
curve slope, and experience has shown that such discontinuities are focus points for further analysis of the stability 
characteristics of the configuration. 
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Figure 16 shows the static rolling moment characteristics of the model at p = 0°. Roll asymmetries are small 
below a = 40°, then build slowly to a large peak between a = 40° and a = 60°. Limited flow visualization studies 
conducted on the transport model suggest that the likely cause of this peak is due to interactions between two vortex 
sheets that roll up over the fuselage and the inboard wings/vertical tail at high a’s. It is noteworthy that the 
wing/vortex interaction appeared to be significantly less than is often observed on fighter configurations - which is a 
possible explanation for the differences in frequency effects previously discussed. Further flow visualization or 
surface pressure measurements are warranted to further investigate these effects. 

Figure 17 shows the classical results from forced oscillation tests conducted at several frequencies and 
amplitudes. The data show no appreciable effects from variations of either frequency or amplitudes up through 
a = 20°. Above a = 20°, a nonlinear variation of the computed derivative occurs, indicating nonlinear roll damping 
characteristics. 

Figures 18-21 show the rolling moment variation as the model oscillates at three different combinations of 
amplitude and frequency while maintaining the same nondimensional maximum roll rate as the model passes 
through <f> = 0° ( p = 0.0188). All four figures are plotted at the same scale to facilitate comparisons. As with the 
fighter model, the data coalesce at a = 0° (Fig. 18) and resemble the ellipses indicative of linear variation of 
Cj and Cjp. At a = 10° (Fig. 19), distortions in the elliptical shape of the data appear, and at a = 12° (Fig. 20) the 

ellipses practically collapse into a single curve, indicative of neutral damping, and exhibit slope changes that are 
amplitude dependent. At a = 40° (Fig. 21) the classical elliptical loops are significantly distorted, with a change in 
direction (twisting) of the loops at the larger amplitudes as indicated by the arrows showing the direction of the 
oscillations. This twisting of the loops at the large oscillation amplitudes is indicative of a change in the damping 
characteristics from unstable at the small amplitudes (highest rates) to stable at the larger amplitudes (lowest rates). 

In general, the non-elliptical shapes of the rolling moment data at higher angles of attack highlight the 
inadequacy of the current linear mathematical models of roll damping. Figure 22 presents the classical forced 
oscillation results for the same data sets as shown in Figs. 18-21. The conventional approach of calculating the roll 
damping parameter involves integrating the out-of-phase (with roll angle) component of rolling moment over the 
entire oscillation cycle. This integration procedure results in an “average” value of the damping parameter over the 
oscillation cycle. The conventional results show unstable damping at the small amplitude oscillation at a=40°, and 
stable damping at the large amplitude oscillation. Inspection of the actual time-history data (Fig. 21) shows all three 
amplitudes of oscillation exhibit unstable damping characteristics at small roll angles. 

The single -point approach predicts a well behaved, stable, and nearly linear variation of rolling moment with 
nondimensional roll rate at a = 0° (Fig. 23). At a = 10° (Fig. 24), the rolling moment variation with roll rate shows 
a slightly nonlinear behavior, and has a slightly lower slope, but is still well behaved and indicates stable roll 
damping. At a = 12° (Fig. 25), the variation of rolling moment with roll rate flattens or becomes slightly unstable. 
At a = 40°, the single-point analysis (Fig. 26) indicates unstable roll damping at small rates and stable damping at 
higher roll rates. This is an indicator of potential wing rock characteristics 7 at this angle of attack; however, more 
data is required at lower rates to fully describe the unstable region - particularly near pb/2V = 0, where a large 
change in rolling moment is seen between the two adjacent points. Observations of the available data at low rates 
seen at the maximum roll angles during the oscillation cycles (Fig. 21) seem to show a stable damping effect. 
Unfortunately, at these roll angles, the combined rotational acceleration and sideslip effects may have an impact on 
the apparent damping characteristics which the single-point approach does not capture since the single -point 
approach includes data only from the mid-point of the oscillation cycle, where rotational acceleration and sideslip 
are zero. 


V. Data Summary and Conclusions 

The data obtained for the low aspect ratio fighter model show strong effects of frequency on the rolling moment 
response to the forced oscillation motions. These results illustrate not only the current limitations of modeling 
dynamic rate effects with the conventional forced oscillation derivative approach, but also of the more detailed 
single point approach, considering that the nonlinear rate effects that resulted from this approach still failed to 
adequately predict the observed characteristics. The strong dependence of the rolling moment characteristics on 
reduced frequency in the stall a range indicates that unsteady and/or nonlinear aerodynamic phenomena are 
important, and that new modeling approaches that account for these effects need to be developed. 
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The conventional damping derivatives calculated from the high aspect ratio transport model dynamic data were 
not adequate to model the characteristics observed during testing. Further analysis showed that the dynamic 
characteristics were nonlinear with roll rate; however, the reduced frequency effects were not evident in the data at 
a’s below 20°. Simple nonlinear mapping of rolling moment coefficient with nondimensional roll rate worked well 
to characterize most of the responses measured in the test. The resulting model represents a vast improvement to the 
current state-of-the art in modern transport simulation and analysis techniques. However, work remains to be done 
in isolating and modeling rotational-acceleration and sideslip effects to better characterize the responses measured in 
the test. 

Both configurations exhibited a requirement to move beyond linear damping derivatives in order to better 
model the aerodynamic effects measured in the roll oscillation tests. The high-aspect-ratio transport configuration 
showed nonlinear roll rate dependence, which is easily adapted into conventional look-up table approaches and can 
be included with only minor changes to simulation architecture. In stark contrast were the results for the low-aspect 
ratio fighter model. The data indicated not only nonlinear dependence on roll rate, but also strong dependence on 
reduced frequency. The variation in aerodynamic response with reduced frequency indicates the possible existence 
of time-dependent nonlinear behavior, which will require additional state variables in the mathematical models. 

Further work should be conducted to identify the causes for different response characteristics between the two 
models. One likely candidate is the wing aspect ratio; however, many other factors such as sweep angle, chord 
length, tail planform and placement among others must also be evaluated. Work in developing math modeling 
techniques to fully describe airplane responses to dynamic motions needs to be pursued. Different motion shapes 
(other than sinusoidal oscillations) and effects of motions about other axes need to be explored as well in order to 
obtain a comprehensive data set over a wide range of conditions. 
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Figure 3. Fighter model mounted in wind tunnel for 
roll-oscillation tests. 



Figure 4. Transport model mounted in wind tunnel 
for roll oscillation tests. 



Figure 5. Characteristic shape for stable static and 
dynamic linear aerodynamics. 
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Figure 6. Static lift data for fighter model mounted 
on roll oscillation test system. 
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Figure 7. Fighter model rolling moment at (3 = 0° on 
roll oscillation test system. 
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Figure 8. Conventional roll damping parameter for 
fighter model from roll oscillation tests at oscillation 
amplitude of 30°. 
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Figure 9. Rolling moment for the fighter model at 
a = 10° during oscillations with p max = 0.04. 
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Figure 10. Rolling moment for the fighter model at 
a = 30° during oscillations with p max = 0.04. 
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Figure 1 1 . Rolling moment for the fighter model at 
a = 50° during oscillations with p max = 0.04. 
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Figure 12. Effect of reduced frequency on rolling 
moment. Fighter model, a = 10°. 
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Figure 13. Effect of reduced frequency on rolling 
moment. Fighter model, a = 30°. 
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Figure 14. Effect of reduced frequency on rolling 
moment. Fighter model, a = 50°. 



Figure 15. Static lift data for transport model 
mounted on roll oscillation test system. 



Figure 16. Transport model rolling moment at P = 0° 
on roll oscillation test system. 
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Figure 17. Conventional roll damping parameter for 
transport model from roll oscillation tests at 
oscillation amplitude of 30°. 
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Figure 18. Rolling moment for the transport model at 
a = 0° during oscillations with p max = 0.0188. 
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Figure 19. Rolling moment for the transport model at 
a = 10° during oscillations with p max = 0.0188. 
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Figure 20. Rolling moment for the transport model at 
a = 12° during oscillations with p max = 0.0188. 



Figure 2 1 . Rolling moment for the transport model at 
a = 40° during oscillations with p max = 0.0188. 



Figure 22. Conventional roll damping parameter for 
the transport model during oscillations with p max = 
0.0188. 
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Figure 23. Single-point results for transport model at 
a = 0° and all frequency-amplitude combinations 
tested. 



-0.15 -0.1 -0.05 0 0.05 0.1 0.15 

pb/2V 


Figure 24. Single-point results for transport model at 
a = 10° and all frequency-amplitude combinations 
tested. 
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Figure 25. Single-point results for transport model at 
a = 12° and all frequency-amplitude combinations 
tested. 
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Figure 26. Single-point results for transport model at 
a = 40° and all frequency-amplitude combinations 
tested. 
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